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Abstract. We argue the behavior of Poynting flux-dominated outflows from 
AGN in the galactic cluster systems by performing three-dimensional MHD sim- 
ulations within the framework of the " magnetic tower" model. Of particular in- 
terests are the structure of MHD waves, the cylindrical radial force balance, the 
(de)collimation, and the stability properties of magnetic tower jets. Transition 
between the jet/lobe and the formation of wiggling jet by growing current-driven 
instability are discussed. 



1. Introduction 



Magnetohydrodynamic (MHD) mechanisms are often invoked to explain the 
launching, acceleration and collimation of jets from Young Stellar Objects, X- 
ray binaries, Acti ve Galactic Nuclei (AGNs), Microquasars, and Quasars (see, 
e.g-M 

eier et al. Il200ll . and references therein). Strongly magnetized jets, par- 
ticularly those with a strong toroidal field encircling the collimated flow, are often 
referred to as "current-carrying" or "Poynting flux-dominated" (PFD) jets. A 
large current flowing parallel to the jet flow is responsible for generating a strong, 
tightly wound helical magnetic field. The global picture of a current-carrying 
jet with a closed current system link ing the magnetosphere of the central engine 
and the hot spots was introduced bv lBenfordl (jl978l . l2006h and applied to AGN 
double radio sources. This closed current system includes a pair of current cir- 
cuits, each containing both a forward electric current path (the jet flow itself, 
with its toroidal magnetic field, toward the lobe), and a return electric current 
path (along some path back to the AGN core). 

Theory of m agnetically driven outfl ows in the electr omagnetic regime has 
been proposed by iBlandford I ( 19761) and Lovelace (Il976h and subsequently ap- 
plied to rotatin g black holes (Biandford Znaiek I Il977l ) and to magnetized 
accretion disks ( Biandford k, Payne 19821 ). An underlying large-scale poloidal 
field for producing the magnetically driven jets is almost universally assumed 
in many theoretical/numerical models. However, the origin and existence of 
such a galactic magnetic field are still poorly understood . In contrast with the 
large-scale field models, Lynden-Bell ( Lvnden-Bell 19961 ) examined the "mag- 
netic tower"; expansion of the local force-free magnetic loops anchored to the 
star and the accretion disk by using the semi-analytic approach. Global magne- 
tostatic solut ions of magneti c towers with external thermal pressure were also 
computed bv lLi et al. 1 {2001) using the Grad-Shafranov equation in axisymme- 
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try. Full three-dim ensional MHD num erical simulations of magnetic towers have 
been performed by iKato et al. 



(2004). 



Recent X-ray and radio observations have revealed the dynamical interac- 
tion between the outbursts driven by AGN and the backgro und IGM/ICM, such 
as X-ray "cavities" with corresponding; radio "bubbles" (iFabian et al. 1 12000b 

" 20071 ). 



Croston et al~l l2003t iKraft et al. 1 120061 : IWise et al. 1 120071 : iGitti et al 



The cavities are believed to be filled with very low density relativistic plasma, 
inflated from the electromagnetic jets that are being powered by AGNs. This 
paper describes nonlinear dynamics of propagating magnetic tower jets in galaxy 
cluster scales (> tens of kpc ) based on three-dimens ional MHD simulations to ar- 
gue the jet/lobe transit ion ( Nakamura et al. 20061 ) and the stability properties 



( Na kamura et al. 1 120071 ) . 



2. Numerical Setup 

We solve the nonlinear system of time-dependent ideal MHD equations numeri- 
cally in a 3-D Cartesian coordinate syst e m (x, y, z). The basic numeri cal treat- 
ments are introduced in lLi et al.~l (|200d ); lNakamura et al. I (|2006l . l2007l ). We as- 



sume an initial hydrostatic equilibriu m in the gra vitationally stratified medium, 
adopting an iso-thermal King model ( King |[l962l ) to model the magnetic towers 



from AGNs in galaxy cluster systems. AGN accretion disk can not be resolved in 
our computational domain, and thus the magnetic flux and the mass are steadily 
injected in a central small volume during a certain time period. Since the in- 
jected magnetic fields are not force-free, they will evolve as a "magnetic tower" 
and interact with the ambient medium. In the present paper, we present two 
different runs: one is called the "unperturbed case" in the following discussion, 
whic h is a run without any initial perturbation to the background initial pro- 
files ( Nakamura et al~ll2006l ). The other is called the "perturbed case", where 



a finite amplitude perturbation (a few percent of the background sound speed) 
is given to the velocities of the background gas ( Nakamura et al."ll2007l ). 



The total computational domain is taken to be \x\, \y\, \z\ < 16 corre- 
sponding to a (160 kpc:— 80 to 80 kpc) 3 box. The numerical grids are 240 3 
in the unperturbed case and 320 3 in the perturbed case. Normalizing fac- 
tors are a length Rq = 5 kpc, a sound speed C s o = 4.6 x 10 7 cm s -1 , a time 
to = 1-0 x 10 7 yr, a density po = 5.0 x 10~ 27 g cm~ 3 . The corresponding unit 
pressure po as poC^ = 1.4 x 10~ n dyn cm -2 , and the unit magnetic field Bq 

as ^AttpqC^q = 17.1 pG. The initial sound speed in the simulation is constant, 

C s = 7 1 / 2 ~ 1.29, throughout the computational domain, which give a sound 
crossing time r s « 0.78, corresponding to a typical time scale Ro/C s o « 10.0 
Myr. Therefore, t = 1 is equivalent to the unit time scale 12.8 Myr. In the King 
model we use here, we adopt the cluster core radius R c to be 4.0 (i.e., 20 kpc) 
and the slope k to be 1.0 in the unperturbed case and 0.75 in the perturbed 
case. Magnetic fluxes and mass are continuously injected into a central volume 
of the computational domain for i; n j = 3.1, after which the injection is turned 
off. A magnetic energy injection rate is ~ 10 43 ergs s , a mass injection rate is 
~ 0.046Mo/yr, and an injection time is ~ 40 Myr. 
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Results 
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Figure 1. Distribution of density in the x — z plane along with the poloidal 
velocity (left) and three-dimensional view of selected magnetic field lines 
(right) for the perturbed case. 



During the dynamical evolution of magnetic tower jet, the narrow "jet 
body" and the expanded "lobe" are formed as shown in Fig. 1 (left). The 
3D view of magnetic field lines as illustrated in Fig. 1 (right) indicates that the 
magnetic tower jet has a well-ordered helical field configuration showing that a 
tightly wound central helix goes up along the central axis and a loosely wound 
helix comes back at the outer edge of the magnetic tower. The profiles of un- 
derlying external gas plays an important role in the transition of jet/lobe. The 
jet body is confined jointly by the external pressure and the gravity inside R c , 
while it expands radially beyond R c to form the lobe in a steeply decreasing 
ambient pressure. 

The jet axial current and the ambien t gas pressure can dete rmine the ra- 



dius of the magnetically dominated lobes ( Nakamura et al. 20061) , that is, the 



internal gas pressure pl ays a minor role in the lobes which are typically seen 
in FR I radio galaxies ( Croston et al. I [20031 ) . Our numerical result shows the 
magnetical ly dominated lobes expand subsonically which can be seen in some ob- 
servations ( Kraft et al. 2006). A preceding hydrodynamic shock wave beyond 
the tower, w hich may be identifi e d as an AGN-drive n shock in recent X-ray 
observations ( Forman et al. |[2005l ; iFabian et al. II2006I ). 

The outer edge of the magnetic tower may be identified as a tangential 
discontinuity without the normal field component. The interior of tower (lobe) is 
separated from the non-magnetized external gas via this discontinuity. Figure 2 
shows the distributions of forces (left) and the magnetic field components (right) 
along the radial directions. At the tower edge, the outwardly directed magnetic 
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Figure 2. Radial profiles of the forces (left) and the magnetic field compo- 
nents (right) along the £-axis for the unperturbed case. The position of the 
magnetic tower edge is shown. 




Figure 3. Distribution of axial current density in the x — z (left) for the 
perturbed case. Time variation of the azimuthal Fourier power of the current 
density with the axial (z) wavelengths A = 6.0 (middle) and A = 2.0 (right). 



pressure gradient force is roughly balanced with the inwardly directed thermal 
pressure gradient force. On the other hand, at the core part of jet body, the 
quasi- force free equilibrium is achieved. In the context of magnetically controlled 
fusion systems, the helical field in the magnetic tower can be regarded as the 
reversed field pinch (RFP) profile. 

The current-carrying magnetic tower jet, which possesses a highly wound 
helical field, is subject to the current-driven instability (CDI). Although the 
destabilizing criteria will be modified by the ambient gas and the RFP con- 
figuration of the tower, we find that the propagating magnetic tower jets can 
develop the non-axisymmetric CDI modes. Both the internal elliptical (m = 2) 
mode (Fig. 3 right) like the "double helix" and the external kink (m = 1) 
mode ( Fig. 3 middle) grow to produce the wiggles at the different parts (Fig. 
3 left) ( Nakamura et al. 1120071 ). This morphological feature (the jet is suddenly 
disrupted at the lobe region with v isible [external] wiggles) can be observed in 
Hercules A (jGizani fe Leahy II2003T ). 
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4. Conclusion 

By performing 3D MHD simulations, we have investigated the nonlinear dynam- 
ics of magnetic tower jets in the galaxy cluster system. Under the gravitationally 
stratified atmosphere, the magnetic tower jet forms a global current closure sys- 
tem having the reverse field pinch magnetic configuration. The magnetic tower 
jet is well confined by the external gas pressure to form the collimated "jet body" 
inside the cluster core radius. The narrow jet expands radially to form the mag- 
netically dominated "lobe" beyond the cluster core radius. Thus, the profiles of 
underlying external gas play a important role in the transition of jets/lobes. 

The current-carrying tower jet, which possesses a highly wound helical mag- 
netic field, is subject to the current-driven instability. In the gravitationally 
stratified atmosphere, the effect of thermal confinement on the magnetic tower 
jet will get gradually weaker as it grows. As a result, a separation of current flow- 
ing path between the jet axial current and the return current occurs. Therefore, 
the edge of the jet axial current becomes a free boundary against the external 
kink mode to produce the visible wiggles as seen in some observations. 
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